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Abstract 

We propose a simple model of classical open system consisting of two subsystems all stationary 
states of which correspond to phase synchronization between the subsystems. The model is gener- 
alized to quantum systems in a finite-dimensional Hilbert space. The analysis of the simplest two 
qubit version of the quantum model shows that all its stationary states are nonseparable. 
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Entanglement of quantum states is a notable resource of quantum information science and 
plays a key role in it's advanced applications such as superdense coding, quantum teleportation 
and quantum cryptography. The important question of whether there is any analogue 
of the entanglement in classical systems remains open. This problem has been considered 
by different approaches 2r|4j, but all of them focused on closed Hamiltonian systems only 
and tried to find a connection between the entanglement of ground state of the system and 
peculiarities in its phase portrait under the variation of parameters. 

The goal of the present paper is to establish a relationship between the entanglement 
of stationary states of an open quantum system and such well-known classical nonlinear 
phenomenon as synchronization. As a first step on this way, we consider a simple model of 
phase synchronization of two subsystems which can be generalized to a quantum case. The 
analysis of the simplest two qubit version of the model leads to the main result of the paper: 
All stationary states of the system turns out to be entangled. Although this result is obtained 
only for a particular model and has not been generally proved, it suggests a new promising 
perspective for better understanding of mixed entangled states properties. 

Consider a simple dynamical model of open system with three variables l x ,l y , and l z whose 
time evolution is governed by the following coupled equations 

~~dt 2l x lyi (1) 

^ = -2U, 
dt 

This system may be considered as the model for phase synchronization of two subsystems 
for the following reasons. First, remind that synchronization in nonlinear dynamics is a 
phenomenon of rhythms adjustment of oscillating systems due to weak interaction between 
them [[I . Here we consider only the simplest case of the phase synchronization when the phase 
difference of two subsystems vanishes with time. Using phases of subsystems (fx, ^2 and their 
amplitudes ri, r 2 as dynamical variables, we can write the complete system of equations in the 
form: dipi/dt = fi(ipi, y?2, r±, r 2 ) and dri/dt = gi(<pi, f2, r i, ^2), where i = 1,2 and /$ and gi are 
nonlinear functions providing the synchronization. Then, this system of four equations can 
be rewritten as a system of two equations by introducing complex variables z\ = r x exp i(fii 
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and z 2 = r 2 exp iip 2 . Finally, with new variables l x , l y , l z connected to zi, z 2 by the relations 



z\z 2 + Z\Z 2 



I, 



i{z l z 2 f - z\z 2 ) 



z l 



Z 2 



(2) 
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we eventually arrive at a system of nonlinear differential equations for l x ,l y , and l z . Note that 
the transformation defined by Eq. (j2j) is a classic analogue of the Schwinger transformation in 
quantum mechanics which used for the description of two Bose oscillators in terms of angular 
momentum components. Since l y can be expressed as l y = r\r 2 sin(<^2 — fi) we assume 
that vanishing of l y is the condition of the phase synchronization between the subsystems. 
Returning to system (TjQ), we note that it has two integrals of motion: I 2 = l x + I 2 + I 2 and 
k = l y /l z . A simple analysis shows that under arbitrary value of I 2 the solution of Eq. (pQ) 



tends to the final state l y — l z — and l x = I. Thus the presence of phase synchronization in 
the classical model represented by Eq. (pQ) is proved. 

Now let us discuss the question about possible quantum analogues of model ([1]). The 
proper way to quantize Eq. (pQ) (at least in semiclassical approximation) has been proposed 
by the author in Ref. [6|. In the present case this method of quantization can be formulated 
as follows. First, we have to represent the classic equations for components l x , l y , l z in the 
form allowing quantization: 



dt 



- lx 



dH 
dl 



+ i x{tR —- tI r- 



(3) 



where H is real and R and R* are complex functions of l x , l y , l z (star means complex conju- 
gation). Then, following Ref. [6], a self-consistent quantum version of Eq. ([3]) can be obtained 
if one writes down the Lindblad equation for density matrix of the system: 



dp 
~dt 



H ,p 



+ 



Rp, i? H 



(4) 



where Ho and R, R + are operator analogues of functions H , R, R* in Eq. <^j. 

Certainly the question regarding the order of operators /, arrangement in H Q and R exists 
because of their noncommutativity, but in semiclassical approximation this problem does not 
arise. It is easy to check by direct verification that system (JTJ can be represented in the form 
required for quantization of Eq. (J3J) by substitution H = f (I 2 ) , R = l z — il y (the first term 
in r.h.s. of Eq. ()3]) vanishes in this case). Thus we see that possible quantum analogues of 
classical model of synchronization ([!]) can be described by the Lindblad equation 



dp 
~dt 



Rp,R A 



+ 



R,pR* 



(5) 



where R = l z — il y . 



It should be noted that using the evolution equation for the density matrix, given Eq. (J5J), 
we can write equations of motion for the average value (A) of an arbitrary observable A. For 
example, if one takes l x as A, the result is 



dll 



dt 



l X) R 



+ c.c. 



l z +ily) [Iz-ily) ) + c.c. 



(6) 



We see that in semiclassical approximation Eq. (JBJ) coincides with the first from Eqs. 
(OQ) as it should be. Our next task is to find stationary solutions of Eq. (jSJ) and to eval- 
uate their entanglement. The simplest situation in which this can be done exactly is a 
two qubit realization of Eq. (j5J). In this case operators l x ,l y ,l z have the following rep- 
resentation: l x = \{a x <g) 1 + 1 <g) a x ), l y = \{(T y ® 1 + 1 <8) a y ) , l z = \{p z ® 1 + 1 <8> a z ), 
where a x ,a y ,a z are ordinary Pauli matrices. Operator l 2 = V x + l 2 y + l\ has the form 
I 2 = | (3 + a x <g) a x + <j y (g) a y + a z (g) cr z ) and commutes with all ^. Operator R = l z — il y 
has the following matrix representation: 

(2 -1 -1 \ 
10 0-1 
10 0-1 
\0 1 1-2/ 



(7) 



With the help of Eq. ([7]) we can find that equation R \ ^/) =0 has two linearly independent 
solutions: 



l*i> 



1 
1 

w 



and 



1^2) 



1 

71 



/o\ 

1 

-1 

V / 



Therefore in the two qubit case the general stationary solution of Eq. ([5]) can be written as 



p st = a + b |* 2 > (*2| + c (^2 1 + h.c. 



(8) 



Coefficients a, b, and c in Eq. ([8]) must satisfy two conditions: a + b = 1 and ab > c 2 which 
correspond to normalization and positivity of matrix p st . Using Eq. ([HD we can write the 
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desired expression for p st as follows 
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It is easy to check that p st has two null eigenvalues and two positive ones, Ai and A2, which 
satisfy to the following equation: A 2 — A + ab — c 2 = 0. Let us find the values of coefficients 
a,b,c, for which density matrix p st (a,b,c) corresponds to entangled states of two qubits. A 
simple way to accomplish this is to invoke the Peres criterium Q. As well known Q], in 
the two qubit case this criterium is a necessary and sufficient condition of the mixed states 
separability. According to it density matrix p is separable if matrix ppx obtained from p by 
operation of partial transposition (which corresponds to permutations of indices of one of the 
subsystems only) is nonnegative. Using Eq. (jSJ) for p st we can write (p st ) PT as 
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It is easy to check that (p s t) PT has eigenvector 

1 

|Ai) 



1 

V2 
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with eigenvalue Ai = b/2. The remaining eigenvalues are found from characteristic equation 



Computing this determinant explicitly we get the following cubic equation for Ai, A 2 , and A 3 

\ 3 -(a + ^\ 2 -(^ + c 2 -f)\ + ^ = 0. (12) 



We see that two roots of Eq. ( TT2T) are positive but the third is negative (when 6^0). Thus 
for all b >- density matrix p st corresponds to nonseparable (entangled) state. This result 
obtained for the simple model of synchronization together with some qualitative reasons 
suggests that a relationship between the phase synchronization in classical open system and 
the entanglement in its quantum analogue exists in more general situations as well. 

In conclusion, we want to examine our original model ([1]) from somewhat different point 
of view. Let us consider two functions of state H = {I 2 . + l 2 + l 2 )/2 and S = 2l x . As follows 
from Eq. (pQ) the evolution of H and 5* satisfies two general conditions: 

dH , dS 

— = and — > 0. 13 

dt dt ~ v 1 

It should be noted that system (pQ) can be represented as: 

dli dH 

Tt = £M dT k Al > (14) 

where A[ = ei mn dS / dl m dH / dl n and ei mn is the antisymmetric tensor of Levi-Civita. 

As shown by the author [9], two conditions (fT3"j) (in the general case of n variables x\, . . . , x n 
describing the evolution of the system) define a class of nonhamiltonian systems called qua- 
sithermodynamic with interesting dynamica. and S tatistica. propert.es g Bntthe possibiHty 
of quantum quasithermodynamic systems existence was not discussed in Ref. [9J. The analysis 
presented above shows, in particular, that a quantum analogue of system ([1]) exists and it is 
a quasithermodynamical system as well. One can see it directly from the fact that average 
values of H = f 2 /2 and S = 2l x obviously satisfy to conditions ( fT3l) . Sure, more detail analysis 
of quasithermodynamic quantum systems is required, which is the subject of future work. 

I thank L.A. Pastur for valuable discussions of results of this paper. 



[1] M.A. Nielsen and I.L. Chuang, Quantum computation and quantum information (Cambridge 

University Press, Cambridge, 2000). 
[2] A.P. Hines, R.H. McKenzie, and G.J. Milburn, Phys. Rev. A 71, 042303 (2005). 



6 



[3] C. Emary, N. Lambert, and T. Brandes, Phys. Rev. A 71, 062302 (2005). 

[4] M.S. Santhanam, V.B. Sheorey, A. Lakshminarayan, Phys. Rev. E 77, 026213 (2008). 

[5] A. Pikovsky, M. Rosenblum, and J. Kurths, Synchronization. An universal concept in nonlinear 

sciences (Cambridge University Press, Cambridge, 2002). 

[6] E.D. Vol, Phys. Rev. A 73, 062113 (2006). 

[7] A. Peres, Phys. Rev. Lett. 77, 1413 (1996). 

[8] M. Horodecki, P. Horodecki, and R. Horodecki, Phys. Lett. A 223, 1 (1996). 

[9] E.D. Vol, physics.class-ph 0812.3738. 



7 



